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We present the results of a NNLO QCD analysis of the World data on unpolarized 
DIS Non-Singlet Structure functions. 



1. Introduction 

The increasing accuracy of DIS experiments will further reduce the experi- 
mental errors on the determination of the strong coupling constant calling 
for an improvement on the theoretical errors, which are by now the dom- 
inant ones. One way to achieve it is to include NNLO QCD effects in the 
analysis. 

Our goal is to perform a NNLO QCD analysis of World data on unpo- 
larized DIS structure functions to determine a s with an accuracy of 0(2%) 
along with a parametrization of the parton distribution functions with fully 
correlated errors. As a first step in this direction we concentrate on the non- 
singlet (NS) sector. 

We presented the first results of our analysis in 1 . Here we give an update 
of the main results and refer the interested reader to 2 for all the details. 



2. Theoretical Framework 

We carried out our analysis in Mellin— N space 3 , where the non-singlet part 
of the electromagnetic DIS structure function F%{N, Q 2 ) is written in terms 
of the non-singlet quark combinations q ,V {N ,Q 2 ) and the corresponding 
Wilson coefficients (N) as 

F^(N,Q 2 ) = {l + a s (Q 2 )C ( 2 1) (N)+a 2 s (Q 2 )ci 2 \N)} q ± ' v (N,Q 2 ), (1) 
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with a s (Q 2 ) = a s (Q 2 )/4:'K, the normalised coupling constant. 

In the region x > 0.3 we adopt the quark valence dominance hypothesis 
under which the proton and deuteron structure functions are given by the 
following quark distribution combinations 

4 1 5 

F% = -xu v + -xd v , F% = — x{u v + d v ) . (2) 

For x < 0.3 we analyse the NS combination 

F 2 NS = 2(F% - F$) = ±x(u v - d v ) - ^x(d - u) . (3) 

The valence parton distribution functions are parametrized at the ref- 
erence scale Qq = 4 GeV 2 with the functional form 

xu v (Ql,x) = A u x au (l-x) K (l + PuVx + Jux) (4) 

and 

xd v (Ql,x) = A d x a «(l-xp(l + Pd ^x~ + ~f d x), (5) 

where the normalization constants A u and A d are not free parameters of the 
fit but are determined to satisfy the valence quark counting: J u v (x)dx = 2 
and L d v (x)dx = 1. 

The remaining non-singlet parton density, (d — u) , is not constrained 
by the electomagnetic structure function data and we adopted the form 
given in 4 , which provides a good description of the Drell-Yan dimuon pro- 
duction data from the E866 experiment. The heavy flavor corrections were 
accounted for as described in 5 . 

3. Data 

The results we present are based on 551 data points for the structure 
function F^x, Q 2 ) measured on proton and deuteron targets. The exper- 
iments contributing to the statistics are: BCDMS 6 , SLAC 7 , NMC 8 , HI 9 
and ZEUS 10 . 

The BCDMS data were recalculated replacing Rqcd with f^iggg 11 . All 
deuteron data were corrected for Fermi motion and off-shell effects 12 . 

We used the measured structure functions F% and F% in the region 
x > 0.3 which is expected to valence dominated, while in the region x < 0.3 
we construct the non-singlet structure function F^ s = 2(F% — F$) from 
proton and deuteron data measured at the same x and Q 2 . 

We imposed different cuts on the data. Only data points with Q 2 > 4 
GeV 2 were included in the analysis and a cut on the hadronic mass of 
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W 2 > 12.5 GcV 2 was imposed in order to reduce higher twist effects on 
the determination of Aqcd and the PDF parameters. The latter cut was 
then relaxed in the extraction of higher twist effects. Moreover we imposed 
additional cuts on BCDMS (y p > 0.3) and NMC (Q 2 > 8 GcV 2 ) data in 
order to exclude regions with potentially significant correlated systematic 
errors. 

In the fitting procedure we allowed for a relative normalization shift 
between the different data sets within the systematic uncertainties quoted 
by the single experiments. These normalization shifts were fitted once and 
then kept fixed. 



4. Results 

The results we obtain for the fit parameters are collected in Table 1. Wc 
note that the fit doesn't constrain the pi and ji parameters, which have 
therefore been kept fixed after the first minimization and their value is 
quoted without errors. 

The remaining parameters to be determined in the fit are the low- and 
high- a; ones (dj and hi) alongside with Aqcd- From the value for ^q CD 



Table 1. Parameters values deter- 
mined in the NNLO QCD fit. 



U v 


a 


0.291 ± 0.008 




b 


4.013 ± 0.037 




P 


6.227 




7 


35.629 


dv 


a 


0.488 ± 0.033 




b 


5.878 ± 0.239 




P 


-3.639 




7 


16.445 




226 ± 25 


X 


z /ndf 


472/546 = 0.86 



obtained in the fit we extract the following value for the strong coupling 
constant 

a s (M 2 ) = 0.1134+°;°°^ (expt). (6) 

We note that this value is in agreement within the errors with results ob- 
tained from other NNLO QCD analyses 13 ' 14 and with the the world average 
0.1182 ± 0.0027 15 within 2a. 

In Figure 1 we compare the parton distribution functions xu v {x) and 
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Figure 1. The parton densities xu v (left) and xd v (right) at the input scale Qq = 4 GeV 2 
compared to results obtained by MRST 13 and Alekhin 14 . The shaded areas represent 
the fully correlated la error bands. 
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